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Under certain conditions, the propagation of a beam of oscillating neutrinos in a variable-density medium takes 
the form of an almost complete transformation of the initial type of neutrino into another type. The depth of 
oscillations is then negligible. The transformation can occur in the cores and envelopes of collapsing stars. 



1. Introduction 

The interaction of neutrinos and matter modifies 
the picture of I'-oscillations P, Si- The effect of 
the medium is analogous in a number of significant 
respects to the coherent regeneration of K mesons, 
and also to the appearance of a refractive index. 

In a medium of constant density, the overall char- 
acter of the oscillations is the same as in vacuum, i.e., 
there is a change in only the oscillation length and 
depth. A number of astrophysical [ll, 0, 0] and geo- 
physical P, d, ^ applications of these oscillations have 
been considered in the constant-density approxima- 
tion. 

In the present paper, we examine qualitatively new 
effects that appear in a variable-density medium. 



2. Equation for the transformation 
probability; resonance condition 

We shall consider the mixing of two types of neu- 
trino. To be specific, we shall suppose that they are 
i/e and Vfj_. Suppose the are created in the source. 
We shall seek the probability P{t) of an oscillatory 

— > transformation in a time t (or distance r ^ ct 
from the source). The equations describing the evo- 
lution of the wave functions of the Ve and in the 
medium [l[ (essentially, the Schroedinger equations) 
can be used to show that the probability P{t) is the 
solution of the following equation: 



A^P dM(PP 



-2M^^{2P- 
dt ^ 



)-- 
' dt 

1) = 0, 



(1) 



where 



M = {2n/l^){cos2e-l^/lo), 2M = (27r/;^) sin 26*, 

(2) 



ll, = AttE/ Am? is the oscillation length in vacuum, E 
is the neutrino energy, Am^ = — rn^ is the differ- 
ence between the squares of the masses, 9 is the mix- 
ing angle, Iq is a characteristic length of the medium, 
given by 



rriN ^ 



A/.(0) 
2-kE 



(3) 



A/j(0) — /je(0) — fifi{0) is the difference between the 
Ve and forward-scattering amplitudes for the ith 
component of the medium (i = e, p, n), Xi is the abun- 
dance of the ith component per nucleon, p is the den- 
sity of the medium, and mjv is the nucleon mass. In 
accordance with (3), we have Iq^ Gpp/'mN, where 
Gi? is the Fermi constant. 

The initial conditions for (1) are: 



P(0) = 1, 



rfP(O) 
dt 



0, 



d^P 
dt^ 



-2AP 



(4) 



If the density is constant (dM/dt = 0), we find from 
(1) that 

This equation has a periodic solution of the form 

P= 1 - Afim^{TTr/lm) 
with oscillation length 



27r (M^ + 4Af 2 



'1/2 



(6) 



and oscillation depth determined by the mixing angle 
in the medium ^: 



A = sin^ 29^ = M^^^/tt^ 



(7) 



^ Apart from dm, the oscillation depth is also found to depend 
on the initial conditions. The formula (7) then corresponds to 
the initial condition (4). 
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We recall that the mixing angle 9m relates the Ve and 
Vij, states with the neutrino eigenstates z^j™ and 
in the medium. (In the medium the states with the 
specific masses and V2 arc not eigenstates of the 
Hamiltonian and themselves oscillate). 

Using (2) and (6), we can rewrite (7) in the form 

sin^ 26l„ = sin^ 20[(cos26' - l^/laf + sin^ 29]-^ (8) 

from which it follows that, for low values of sin^ 29, the 
dependence of sin^2^TO on lu/lo{lu/lo ~ pE) exhibits 
a resonance. When 

1^/Iq = cos29, (9) 

sin^2^TO is a maximum, i.e., sin^ 2^^ = 1. The 
quantity sin^ 20„i falls rapidly as l,y/lo departs from 
cos2(^ = 1. Condition (9), for which the mixing angle 
in the medium is equal to 45° (sin^ 29m = 1), will be 
referred to as the resonance condition. Correspond- 
ingly, the values of p and E for which (9) is satisfied 
will be referred to as the resonance values. From (9) 
and (3), we have 

PR = mN^m'^cos2e/2^/'^GFE. (10) 

The width Apu of the resonance layer will be defined 
as the density interval around pR in which sin^ 29m > 
1/2. Prom (8) we have 

Apn = pR,tSin29. (11) 

The resonance becomes narrower as the mixing angle 
decreases. Similarly, we may introduce a resonance 
energy Er and a resonance width AEr = Er tan29. 

The physical meaning of the resonant behavior of 
sin^ 29m is as follows. Suppose that a constant-density 
layer intercepts a neutrino flux with a continuous en- 
ergy spectrum. Neutrinos with energy E = Er{p) 
will then oscillate with maximum oscillation depth 
(despite the small mixing in vaccum). The energy 
dependence of the oscillation amplitude is determined 
by the resonance curve (8): sin^ 29m = .f ['i^ (£")]• 

In the case of a variable-density medium, the oscil- 
lation depth does not depend on 9m alone, and even 
in the resonance layer (p = pr), it may turn out to be 
low (see below). 

3. Medium with slowly-varying density 
(adiabatic regime); nonoscillatory 
transformation in the u-beam 

Consider a layer of a medium satisfying the follow- 
ing conditions: 

(1) The density distribution p{r) has no singulari- 
ties and the derivative dp/dr is a smooth function. 

(2) The variation of density with r is relatively slow 
and such that, for a given mixing angle, the width Atr 



of the resonance layer in the r scale is larger than the 
oscillation length 

Atr = {dp/dr)-^ApR > 1^/2. (12) 

This will be referred to as the adiabatic condition. 

(3) The density distribution p{r) is a monotonic 
function of r up to pmin ~ 0. We shall also assume 
that the resonance density for neutrinos of given en- 
orgy E falls into the interval po — Pmin, so that the 
neutrinos generated in the region with p = po will 
traverse the resonance layer. 

Under these conditions, the solution P{r) is uni- 
versal with respect to the density distribution. The 
universality can conveniently be expressed in terms of 
the following dimensionless parameters. Let 

n = [p{r) - Pr]/Apr, (13) 

where pr and ApR are determined by the values of E, 

9, and Am? . We shall measure distance in units of n 
rather than r. We note that, at resonance, n = 0, and, 
as /9 ^ oo, n ^ oo, while for p — > 0,n ^ — (tan2^)~^. 
The initial conditions are set at 

no = [p(ro) - Pr]/Apr, (14) 

where hq is the number of resonance layers that can 
be fitted between the point at which the neutrino is 
created and the resonance layer. When p(ro) '> Pr, 
we have 

no « p{ro)/pRsm29 w 1/ sin 26*^ 

[see also (8)], i.e., no is equal to the reciprocal of the 
mixing parameter at the point at which the v is cre- 
ated, and increases with distance from the resonance 
layer and/or with decreasing sin^ 29. 
We also introduce the further variable 

m = {dp/dr)-' Apr/ (15) 

i.e., the number of oscillation lengths in the resonance 
layer. This number increases as the distribution p{r) 
becomes shallower or sin^ 29 increases. 

When conditions (l)-(3) are satisfied, the solu- 
tion P{r) is a function of n, uq and m, i.e., P ~ 
P(n, no, m), and is not very dependent on the density 
distribution p{r). 

Let us examine some of the properties of the func- 
tion P{n, no, m): 

(a) P is an oscillating function of n with period T; 
by definition, T ^ 1/m in the resonance layer and, as 
we depart from resonance, T ~ l/[m(n^ -|- 1)"^^^]. 

(b) P oscillates around its mean value 

P(n,no) = [l + no(n2 + l)-V2„(„2 + i)-i/2]/2; (16) 

(c) The amplitude of the P oscillations is Ap = 
\Pmax — Pmin\ and is a maximum in the resonance 
layer 

A««(ng + 1)-V2 (17) 



Sov. Phys. JETP 64 (1) July 1986 



and decreases with distance from the layer. Ap is 
practically independent of m. 

(d) At exit from the layer, for p ^ 0, 

Po-^[l-m{nl + l)-^/^cos2e]/2. (18) 

Consider the limit of large no- Increasing no means 
that the point at which the v is created becomes more 
distant (in p) from the resonance layer. The oscil- 
lation amplitude is then shown by (17) to decrease 
{Ap ~ I/uq) and the average of P(n, no) tends to the 
asymptotic value given by 

P„(n) = [l + n(n2 + l)-V2]/2. (19) 
At exit, Pa sin^ 6 (see Fig. 1). 




Figure 1: Oscillation probability P{n) for different 

values of no (indicated against the curve) . Solid curve - 
asymptotic form of P{n) (nonoscillatory transformation). 

Thus, for small mixing angles in vacuum and 
large no (so that the v-flux. is generated quite far from 
the resonance layer), the neutrino propagation process 
takes the form of a virtually nonoscillatory transfor- 
mation of one type of neutrino into the other. 

The interpretation of these results is as follows. The 
neutrino oscillations in the medium take place around 
the cigcnstatcs z/™ and , i.e., the neutrino oscillates, 
as in the vacuum, around states of specific mass of vi 
and 1^2- If the density of the medium varies, there is 
also a variation in the eigenstates z^™ and f ™ or, more 
precisely, there is a change in the mixture of and 
Vn- When p p^, it can be shown from (8) that 
Om ~ 7r/2, i.e., i/™ becomes practically identical with 
fg. When p = pu, wc have dm = vr/4 and i/^ contains 
equal admixtures of and t'p. When p <ti pr, vfe have 
9m = and, when 9 is small, v'!^ consists mostly of Vfj,. 
When the density changes from p ^ pji to p <^ pa, 
the basis v]^ , rotates through 90° relative to the 
basis i/f^. The adiabatic condition then shows that 
the neutrino state z^(t) is altered as a result of the 
variation in density: ^{t) follows i/J" and f™, while 
the admixtures i^j" and in u{t) undergo only a 
small change. 



Let us suppose that the neutrinos Ve are gener- 
ated in a region with po ^ Pr and then cross layers 
with continuously decreasing density {p 0). Ini- 
tially, 6l,"„ - 90°, sin^ 29^^ - sin^ 20(pr/ paf < 1, and 
z/(0) = f g ~ i/^. As the density tends to zero, the 
state rotates through 90°, as noted above, and be- 
cause of the adiabatic property, i^™ rotates together 
with v{t). When 9 is small, v{t) ~ vl^ ~ in the 
final state. Thus, becomes transformed into v^. As 
the distance of the point of creation of the neutrinos 
from the resonance layer increases (in p), sin^ 20^ de- 
creases, and the difference between v{t) and j^™ also 
decreases. There is an attendant reduction in the os- 
cillation depth. In the limit of large po, the function 
v{f?j becomes practically identical with z/™, i.e., the 
eigenstate in the medium, and, consequently, there 
are no oscillations. This limiting case corresponds to 
a nonoscillatory transformation. 



4. Applications of oscillation effects 

We shall now formulate the general conditions for 
the above effects to produce observable consequences. 

(1) Both amplification of oscillations and signifi- 
cant changes in the properties of the v-heam. in the 
adiabatic regime arc due to the crossing of the res- 
onance layer by the neutrinos. The resonance con- 
dition is satisfied for a particular sign of Iv/h or 
A/(0) cos 2f?/Am2 [scc(8)]. Since replacing ;/ with v 
produces a change in the sign of A/(0), the resonance 
condition is satisfied in a given medium only for a neu- 
trino or an antineutrino. If the oscillation effects are 
amplified in the i/-channel, they are suppressed in the 
i>-channel, and vice versa. 

(2) the matter-effect occurs in the charge- 
asymmetric medium. If the particle and antiparticle 
densities arc equal, we have SiniA/i(0) = 0. 

(3) The typical scale over which the influence of 
matter is significant is / > ^o = A/Gpp, where A = 
3.5 X 10'* km. Hence, it follows that the thickness of 
the medium must be greater than d = Ip > 3.5 x 
lO^g/cm^. 

These conditions and the adiabatic condition are 
satisfied in the sun as well as in the envelopes and 
cores of collapsing stars. 



5. Neutrino fluxes from collapsing stars 

The effects examined above can occur in the wide 
range of values of Am? and sin^ 26 in the outer layers 
of the cores (above the neutrino sphere) and in the 
envelopes of collapsing stars. 

For the purposes of estimates, we shall use the den- 
sity distribution in the envelope at the beginning of 
the collapse, by analogy with the situation prevailing 
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in white dwarfs: 

P«Po[(i^B/r■)-l]^ (20) 

where p ~ (8 - 10) x 10^ g/cm^ and Rb =5x 10^ cm. 
The maximum density in the interior is Pmax — 10^ 
g/cm^. For massive stars, there is, in addition to (20), 
an extended hydrogen envelope with p « 10~^ g/cm^. 
The envelope may begin to expand during the collapse 
process with velocity v ~ 5000 km/s. 

For collapsing cores, we shall use the model density 
distribution and the neutron and electron densities 
given in Ref. 7. The density in the neutrino sphere 
will be taken to bo p ~ 10^^ - 10^^ g/cm^. 

We have used these density distribution to calculate 
the z^-parameter ranges (see Fig. 2) for which the adi- 
abatic transformation conditions (1) - (3) of Section 
3 are satisfied. The upper limits for Am^ are deter- 
mined by the maximum density. The lower limits for 
sin^ 26 follow from the adiabatic conditions. The val- 
ues of Am^^j, for ^ i/^ and Ug {vg is the 
sterile state) are different: 

A"^ma^(^^e ^ J^^) = 3 • lO'^ Am^„^ (j/g ^ l^s)- 

This difference is largely due to the strong neutroniza- 
tion of the core (n„ He). We note that in the case of 
the ^ i^s)-oscillations, the ratio I^/Iq changes sign 
for p ~ 10'' g/cm^ in the core. Here, we have a cancel- 
lation of the effects on neutrons, on the one hand, and 
on electrons and protons, on the other. This means 
that, if most of the conversion in the core is for the 
neutrinos, most of the conversion in the envelop will 
be for the antineutrinos. 

Since, under the above conditions, the neutrino 
state is practically the same as one of the cigenstates 
in the medium, there is no loss of coherence due to 
the spreading of the wave packets corresponding to 
the different i^™. 

For the ranges of values of the parameters 
Am^,sin^20 given in Section 3 (Fig. 2), the factor 
representing the suppression of the flux of neutrinos 
of the original type approaches something between 
(1/4) sin^ 261 and (1/2) sin^ 261. 

We must now consider the consequences of the os- 
cillation effect in matter from the point of view of 
detection of //-bursts due to collapses. 

(a) Suppose that conditions (l)-(3) are satisfied in 
the z>e f'a channel, where Pa = i'fi or Pt-. In that 
case, the flux will be almost completely transformed 
into the flux, and vice versa; and Da will ex- 
change; their spectra. It is assumed that the 9^ and 
Pt have harder spectra than the Df. (the total fluxes are 
roughly equal for the entire burst). The exchange of 
the spectra between the z>,. and /y„ means that there is 
a substantial increase in the number of events in scin- 
tillation counters detecting Peif^v ~ E^) (Ref. [8]). 

(b) If the resonance condition is satisfied for the 
neutrinos, Va {i^a = ^'/ld^'t)) an effect analogous 




... ■!_ . _J I i \ i I ■ 



Figure 2: Range of (/-parameters in which the cores and 
envelopes of collapsing stars exhibit strong adiabatic 
transformation: a), i/e <-» fp, b). z/e <-» Vs oscillations. 
Below the shaded line-effect at exit from the envelope: 
Po((sin^ ^)/2, dotted line - effect of an expanding 
envelope at the end of the z/-burst t = 10s); dot-dash 
lines - experimental limits. 

to that just described may be expected in systems 
based on radiochemical methods {ve detection). 

The Va i^e or Da i^e spectrum exchange in the 
stellar core leads to a substantial release of energy in 
the envelope because of i/gfi— or D^e— scattering. This 
may play a definite part in the mechanism responsible 
for the shedding of the envelope. 

(c) Resonance in the z/g Vs or D^, Ds channels. 
The sterile neutrino fiuxes appear to be very low. The 
effect due to the precession of spin in the magnetic 
held is negligible for H < lO^^ G and Ato^ < lO^eV^ 
The Vg >^s or D^ — s- Ug transformations will therefore 
give rise to strong (by several orders of magnitude) 
suppression of the v,, or Dg fluxes. 

We emphasize the importance of simultaneous ex- 
periments on the detection of and z/,.. 

We note that, since Iv/Iq has opposite signs in cores 
and envelopes, a strong eff'ect may be present for both 
Df. and z/g with different energies. 

In the case of mixing of three or more neutrino 
types, it is possible that two or more differences Am^ 
will fall into the range of strong adiabatic conversion. 
A combination of the above effects may then be ob- 
served. 

We note that matter effects may not remain 
constant in the course of a v- burst At 5 — 20 s). 
This time dependence will reflect the variation in the 
structure of the core and envelope. 
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Translated by S. Chomet 



Comments (June 2007) 

1. This paper presents, in particular, our first ana- 
lytic results on the adiabatic conversion of neutrinos 
in matter. It has been written in summer-fall 1985. In 
attempt to avoid problems with publication (we had 
before) , we tried to hide the term "resonance" , and 
did not discussed applications to the solar neutrinos; 
also we have not included references to our previous 
papers on the resonance enhancement of neutrino os- 
cillations. 

This short paper has been submitted to JETP 
Letters in the fall 1985 and successfully ... rejected. 
It was resubmitted to JETP in December of 1985. 
The results of the paper have been reported at the 
6th Moriond workshop in January 1986 and included 
in several later reviews. The paper was reprinted in 
"Solar Neutrinos: The first Thirty Years", Ed. J. N. 
Bahcall, et al., Addison- Wesley 1995. 

2. The differential equation of third order for the 
survival probability P, Eq. (1), has been derived 
from a system of three differential equations for P, 
R = Re{i'e\i^fj,) , and I = Im{ye\Vfj). The system of 
equations has been obtained in our first paper: Sov. 
J. Nucl. Phys. 42, 913 (1985), (Yad. Fiz. 42, 1441 
(1985)). 



3. Analytic results of sec. 3, have been derived 
neglecting the high order derivatives d^P/dt^ and 
d'^P/dt^ in Eq. (1) which is implied by the adiabatic 
condition. The resulting equation, 

M(M2 + 4M2)^ - 2M2^(2P - 1) = 0, 



dt 

can be easily integrated: 



P=- 
2 



Pn - - 



dt 



no 



With the initial condition 



1 



P(0) = l--sin^2C = l- 



2{nl + 1) 



it leads to the adiabatic conversion formula (16). 
Noticing that 



cos 26 n 



no 



cos26li 



one realizes immediately that Eqs. (16) and (18) coin- 
cide with the adiabatic formulas that usually appear 
in literature. 
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